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We review and discuss the role of diffeomorphism symmetry in quantum gravity models. Such 
models often involve a discretization of the space-time manifold as a regularization method. Generi- 
cally this leads to a breaking of the symmetries to approximate ones, however there are incidences in 
which the symmetries are exactly preserved. Both kind of symmetries have to be taken into account 
in covariant and canonical theories in order to ensure the correct continuum limit. We will sketch 
how to identify exact and approximate symmetries in the action and how to define a corresponding 
canonical theory in which such symmetries are reflected as exact and approximate constraints. 
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I. INTRODUCTION 

The dynamics of general relativity is deeply inter- 
twined with space-time diffeomorphisms. This becomes 
apparent in the canonical formulation where Hamilto- 
nian constraints generate the dynamics and at the same 
time space-time diffeomorphisms 1 . Moreover requiring 
a canonical theory of spatial geometries and a represen- 
tation of the Dirac algbra - the canonical incarnation of 
diffeomorphisms - results uniquely^ in canonical general 
relativity. 

Hence a quantum theory of geometries with a quantum 
representation of diffcomorpisms should yield a quantum 
theory of gravity with the correct semi-classical limit. 
Unfortunately such a representation seems to be out of 
reach yet. One reason out of many is that many models, 
such as Regge calculus^ and spin foam models^, involve 
an explicit discretization of the space-time manifold as a 
regularization method. Generically this seems to result in 
a breaking of the continuum diffeomorphism symmetries. 

This can be taken as an advantage as such symmetries 
lead to divergencies in the path integral and constraints 
in the canonical formalism and complicate therefore the 
construction of a quantum theory. The diffeomorphism 
symmetries can however not be completely ignored as 
these have to be recovered in the continuum limit. Even 
in discrete theories there should exist at least approx- 
imate symmetries, which should become exact in the 
limit. 

Path integrals for constrained theories such as general 
relativity can be understood as providing a projector on 
the space of physical, that is gauge invariant, states of 
the canonical theory. These are states that are annihi- 
lated by the constraints. The heuristic idea is that path 
integrals implement averagings of the wave functions over 
the gauge orbits. This can only be made consistent if the 
symmetries of the canonical and the covariant theory can 
be brought to match. In particular approximate symme- 
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tries in the action should correspond to some notion of 
approximate constraints. 

We therefore start in section [IT] with a review and dis- 
cussion of symmetries for discretized actions, in particu- 
lar the Regge action. We will clarify that it is not suffi- 
cient to identify infinitesimal variations leaving the action 
invariant. A further requirement is that the associated 
transformations have to act non-trivially on the extrema 
of the action. We will find that in general the symmetries 
defined in this way are a mixture of exact and approxi- 
mate ones. In section IIIII we discuss additional symme- 
tries and problems that can arise if one includes a sum 
over triangulations into the path integral. 

In section IIVI we will line out how to obtain canoni- 
cal formulations that match the symmetries of the dis- 
cretized covariant theory. To this end we use the 
method of 'consistent discretizations' 6 together with a 
local evolution scheme for the dynamics of triangulated 
manifolds 7 . In this way we can make the idea of approxi- 
mate constraints precise. As a simplified example we will 
consider 3d Regge calculus with and without cosmologi- 
cal constant in section IIV Al 

Section [V] discusses notions of Dirac's algebra for dis- 
crete manifolds. In particular we point out that exact 
versions exist for 3d gravity and for 4d theories describ- 
ing flat space. Since 'chunks of flat space' can also ap- 
pear in the full 4d case and moreover in the continuum 
limit curvature should become small compared to the dis- 
cretization scale such a finding could ensure that an exact 
4d Dirac algebra can be recovered in the continuum limit. 

In section IVII we review the role of diffeomorphsim 
symmetries in loop quantum gravity and spin foam mod- 
els and sketch consequences of our earlier discussions for 
these theories. We conclude with a discussion and out- 
look in section |VlTJ 



A. Regge calculus 

The examples we will consider are taken from Regge 
calculus. We will therefore give a very short introduction 
and fix notations. More details can be found in^ and 
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references therein. 

Regge calculus is usually considered on a fixed trian- 
gulation T of a space-time manifold. The variables ap- 
pearing in the Regge action are given by the edge lengths 
{( e } e£ r- This completely specifies the geometry of the 
triangulation. From the edge length one can compute 
the 4d dihedral angles Q° and 3d dihedral angles For 
a 4-simplex a the dihedral angle 0° gives the (inner) an- 
gle between the two tetrahedra sharing the triangle t. A 
similar definition applies to the 3d angles in a tetrahedron 
t at an edge e. 

Consider a triangle t in a 4d triangulation. This tri- 
angle is shared by several 4-simplices. A (Levi-Civita) 
parallel transport of a vector from one 4-simplex to the 
next around the triangle results in a rotation of this vec- 
tor by the so-called deficit angle e* = 2tt — J2a3t ®t m 
the plane perpendicular to this triangle. The deficit an- 
gle measures the curvature concentrated at the triangle. 
Accordingly the (Euclidean) Regge action with cosmo- 
logical constant A is given by 

S = - J2 A t e t + \J2 V ° W 

tebulk a 

where A t denotes the area of a triangle t and V a the 
volume of a 4-simplex a. We work in units with c = 
8itG Newton = 1- If there is a non- vanishing boundary 
the boundary term 

S bd ry = - ]T ( 2 ) 
tgbdry <x3i 

has to be added to the action in order to make the bound- 
ary value problem (with prescribed edge lengths) well de- 
fined. The equations of motion can be obtained by vary- 
ing the Regge action with respect to the length variables 
and are given by 

-i£ et / e cota*+A£^=0 (3) 

where a' is the angle in the triangle t opposite the edge e. 
Derivatives of the deficit angles vanish due to the Schlafli 
identity ensuring that J^tea At89f = for variations of 
the dihedral angles in a 4-simplex. This is analogous to 
the continuum where for the variation of the Einstein- 
Hilbert action J Jgg ah Il a } ) $'x the term with the varia- 
tion of the Ricci tensor leads to a total divergence. 
For 3d Regge calculus we have the action 

S = - he e + Xj2Vr+ £ l e (n-J2 e) (4) 

eGbulk r e£bdry r3e 

yielding the equations of motion 

Here e e = 2n — J2tbc ®l is the 3d deficit angle. If the cos- 
mological constant A = the equations of motion require 
the deficit angles to vanish, hence in this case solutions 
are given by flat space triangulations. 



II. BIANCHI IDENTITIES AND SYMMETRIES 
OF THE REGGE ACTION 

Regge calculus seems to circumvent the problem of 
choosing coordinates and to deal directly with discretized 
geometries. Hence the title of Regge's paper— 'General 
relativity without coordinates'. This however does not 
guarantee that gauge degrees of freedom will not ap- 
pear. Also in the continuum the problem is not so much 
in choosing coordinates, but that different metric fields 
which seem naively to be expressed in the same coordi- 
nates describe the same geometry. 

That gauge degrees of freedom appear in Regge calcu- 
lus was first shown by Rocek and Williams 8 . In a per- 
turbative expansion of the 4d path integral around flat 
space non-propagating modes have been identified and 
it was realized that these correspond exactly to a dis- 
cretization of the continuous gauge degrees of freedom. 
These modes have a simple explanation: any triangula- 
tion of flat space is a solution of Regge's equation. An 
arbitrary vertex of such a triangulation can be translated 
in four directions, keeping the triangulation flat but in- 
ducing the length of the adjacent edges to change. These 
are exactly the gauge modes. 

For a general non-flat triangulation the issue is nat- 
urally more involved. Many author s 8 ' 9 ' 10 ! 11 pointed out 
that there should be at least an approximate notion of 
symmetry if the deficit angles (determined by the edge 
length compared to the curvature radius) are small. The 
arguments rest on the fact that the (linearized) Bianchi 
identity holds approximately, which can be used to show 
that the Regge action should be approximately invariant 
under infinitesimal vertex displacements. 

Some authors argue even for the existence of exact 
symmetries of the Regge action, notably Hamber and 
Williams^. 

Instead of trying to define what the equivalent of a 
continuous diffcomorphism might be for a discrete mani- 
fold, we will take the Regge action and hence the discrete 
system it describes as fundamental and just ask for sym- 
metries of this action, as is advocated for instance by 
Williams 13 . That is we look for infinitesimal transforma- 
tion of the length variables that leave the Regge action 
invariant. 

Indeed one can convince oneself by a counting ar- 
gument that for n-dimensional triangulations at least 
n infinitesimal such transformations per vertex should 
exist. (A similar argument applied to 3-dimensional tri- 
angulations appeared in 8 -.) The piecewise linear man- 
ifold conditions for an n-dimensional triangulation re- 
quire that every point has a neighborhood homeomor- 
phic to the open ball in R™. From this it follows that a 
vertex in such a triangulation is shared by at least (n+l) 
edges. Hence there should be at least n variations of the 
associated length variables that leave the Regge action 
invariant. 

For the 3-dimensional Regge action without cosmolog- 
ical constant such variations were found by Freidel and 
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LouapreJ* with the help of the exact Bianchi identities. 
These transformations were actually defined in a first or- 
der formulation, but can be adopted to 3d Regge calculus. 
We will repeat the argument here, since it will shed some 
light on the 4d case. 

To this end consider the simplest case, that is a vertex v 
that is shared by 4 edges e,, i = 1, . . . , 4. We will choose 
some Cartesian coordinate system, which we can asso- 
ciate for instance to the tetrahedron T%, that has e\,e% y €3 
as its edges. The deficit angles are encoded in parallel 
transport matrices around the edges U\,...U± € 50(3). 
Here U\,U2,Uz start all in the tetrahedron T\ and are 
expressed in the Cartesian coordinates associated to this 
tetrahedron. With U4 we denote the parallel transport 
matrix around e4 starting and ending in some other tetra- 
hedron but conjugated with the parallel transport to T\. 
In this way all parallel transports can be made to start 
and end at the same point and are expressed in the same 
coordinate system. We will denote the resulting matrices 
by Uf . As one can check easily the composition of the 
loops around the edges e\ to can be contracted to a 
point without crossing any edges. Hence the associated 
parallel transport matrix must be trivial, that is we have 

HUi= J]exp(Ai) = l (6) 

i i 

where we defined the Lie algebra elements Aj := \ogUi 
(assuming that the deficit angles arc so small that the 
Ui are in a neighborhood of unity where the exponential 
function is invertible). We will need the Bianchi identity 
([51) in its logarithmic form, however we have to face the 
problem that 3d rotations do not commute. This is the 
point where for instance Miller et al. 10 apply the 'Abelian 
approximation' for small deficit angles 

E A ^° - ( ? ) 

i 

and hence derive an approximate invariance of the ac- 
tion. Freidel and Louapre note that one can use the 
Baker-Campell-Hausdorff formula in order to find the 
exact logarithmic expression for Q in the form 

^(^ + [a,A l ]) = (8) 

i 

where fii are Lie algebra elements depending in a com- 
plicated way on the elements Aj and therefore the deficit 
angles e^. Representing so(3) Lie algebra elements as 
B ab = B c e cah with e cab the totally antisymmetric Levi- 
Civita tensor, ([5]) turns into 

Y^iA' + Afnie^) . (9) 

i 

To find the associated invariance of the action we ex- 
press the relevant part of the Regge action with the help 
of the Cartesian coordinates introduced before as 

s^Y, e - b - A f< ( 10 ) 



where e\ are the edges expressed in the Cartesian coor- 
dinates and Af b — eie^ej/k. (For the edge we use 
the coordinates in T±.) Variations b~e\ of the coordinates 
e\ will induce variations in the associated lengths and we 
obtain 

SS cx £ abc Af 5e1 . (11) 

i 

The terms with SAf = 5(ei)e abc e?/li + e i e abc S(e c l /l i ) van- 
ish. The first summand is orthogonal to ef and for the 
second summand one can use the Schlafli identity. 

With the approximate form of the Bianchi identity we 
obtain that under the variations 

fiapprox^i • — % (1^) 

describing the displacement of the vertex v, the action is 
approximately invariant 

SapproxS = e abc x c S A® b ~ . (13) 

i 

But under the variations 

Sexacte c z := x c - nfx d (14) 
the action is exactly invariant: 

SexactS CX ^ t^Af ( X C - Slf X d ) 

i 

oc J2 2xd ( A t + A i Q t eCed ) 

i 

- (15) 

due to the exact Bianchi identity ([9]). Note that the 
'exact' variations depend via the Lie algebra elements fl; 
on the deficit angles and therefore on the length variables 
whereas the 'approximate' variations do not. In the case 
of a flat vertex, where the = the two transformations 
agree. 

As the counting argument presented above suggest it 
should also be possible to find variations that leave the 
3d action with cosmological constant term exactly invari- 
ant (or any other action for instance with higher order 
terms). For the 4d Regge action the exact (contracted) 
Bianchi identities could be used in a similar way as in 
3d to derive invariant variations corresponding to vertex 
displacements, for alternative actions again the count- 
ing argument applies. Indeed invariant variations can be 
always found. 

The problem is however that the transformations de- 
fined by these invariant variations might turn out to act 
trivially on solutions of the Regge equations, where by 
definition all variations vanish. Consider the example of 
an action depending on two variables and invariant under 
2d rotations in these variables. Due to rotational invari- 
ance there will be an extremum at the center of rotation, 
hence this solution is a fixed point of the symmetry of the 
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action. There could be of course other additional sets of 
solutions with a non-trivial action of the symmetry. 

To test whether candidate symmetries act non- 
trivially on solutions one can consider the Hesse matrix of 
second derivatives of the Regge action evaluated on these 
solutions possesses vanishing eigenvalues or not. The 
Hesse matrix also appears in the perturbative treatment 
of quantum Regge calculus. For the expansion around 
flat space solutions one can indeed find zero eigenmodes 
corresponding to vertex displacements ^ 15 ' 16 . 

To study an (accessible) example with non-flat so- 
lutions, we consider the 3d Regge action with non- 
vanishing cosmological constant. With the additional 
volume term the exact invariance (|14[) does not need to 
hold any more, but as discussed before there should ex- 
ist three independent variations per vertex that leave the 
action exactly invariant. In the case that these vari- 
ations lead to non-trivial transformations also for solu- 
tions, we should find zero eigenvalues for the Hesse ma- 
trix evaluated on these solutions. To simplify the prob- 
lem as much as possible we consider the boundary of an 
equilateral tetrahedron with edge length a and subdi- 
vide this tetrahedron by connecting one inner vertex v to 
the four vertices of the tetrahedron by edges with length 
li,i = 1, . . . , 4. We have to find the extrema of the action 
as a function of the U. There exist one extremum where 
all the inner edge lengths are equal U — b. The length b 
can be found by solving 



, n n 2b 2 -a\ , 

(2tt - 3 arccos — 5 5- J + A 

46 2 - a 2 



a 2 b 



4V36 2 - a 2 



(16) 



For A = we obtain b = y / 3/8a. The equation cannot 
be solved exactly for non-vanishing A but one can ob- 
tain a perturbative solution in powers of A. The (4 x 4) 
Hesse matrix around a configuration with ij = b can be 
calculated 16 and has a very simple form due to the sym- 
metry of our problem. All the diagonal entries are equal 
to a rational function 
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V36 2 - a 2 



2b 2 - 
Ub 2 - 



A 12k 4 - a 4 
12 36 2 - a 2 



(17) 



Also all the non-diagonal elements coincide and are given 
by 



B = 



V3b 2 - a 2 



2b 2 



+ A 



b 2 2b 2 



Ab 2 -a 2 4 3& 2 - a 2 



(18) 



The Hessian matrix will have three null eigenvalues if 
and only if A = B. 

With the equation of motion (|16|) and the requirement 
that A — B we have two conditions for the variable b, 
which can only be satisfied simultaneously for vanishing 
cosmological constant A = 0. For non-vanishing A the 
Hesse matrix evaluated at a solution of the equation of 
motion will not have zero eigenvalues, from which we 
conclude that the candidate symmetries of the action act 
trivially on (at least these) solutions. 



In 4d (without cosmological constant) there are zero 
modes of the Hessian at least for every flat vertex. To the 
best knowledge of the author there is no explicit test for 
vertices with curvature in the A = case, but numerical 
results 17 (involving convergence properties of constraints 
for flat and non-flat solutions) indicate that generically 
there are no zero modes. 

Then different solutions might have symmetry orbits of 
different size, according to whether they have flat vertices 
or not. As will be explained in the next section, there are 
many curved solutions with such flat vertices, moreover 
the number of such flat vertices might become large for 
very fine grained triangulations. 

Additionally one would expect that with fine graining 
and hence decreasing deficit angles gauge symmetry will 
be restored also for the non-flat vertices. In our example 
the only scale available is the cosmological constant A 
(with A = the Regge equations are invariant under 
global rescalings) hence this corresponds to the A — * 
limit. The Hessian matrix will exhibit a set of smaller 
and smaller eigenvalues. The corresponding eigenvectors 
generate an approximate gauge orbit, along which the 
action is almost constant. 

Note that the inverse of the Hessian appears in sad- 
dle point approximations to the path integral. Indeed 
in the same way diffcomorphism symmetry is gradually 
restored the path integral should diverge due to the in- 
tegration over the approximate gauge orbits. Along such 
an orbit the true extremum (or even several extrema if 
the solution is not unique) is difficult to locate. Even 
for classical numerical treatments solving for the exact 
extrema of the action becomes so complicated, that in- 
stead one (gauge) fixes four length variables per vertex 
and does not take the corresponding number of equations 
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into accoun 



III. VARYING THE TRIANGULATION 

One way to quantize Regge calculus is to fix a triangu- 
lation and to perform the path integral over all allowed 
length assignments. One criticism towards this procedure 
is that the final results may depend on the choice of tri- 
angulation, which is rather an auxiliary structure. If the 
triangulation serves as regulator for the continuum path 
integral the induced continuum path integral measure 
might be concentrated around geometries which arise by 
assigning roughly the same length to all the edges 1 ^. 

The dependence on one choice of triangulation can be 
avoided by summing over all triangulations. That is one 
considers a configuration space C labeled by triangula- 
tions and length assignments to the edges of the trian- 
gulations (or more generally some set of variables on the 
different k-simplices) . Also the space of triangulations 
needs to be specified in more detail, for instance whether 
one considers fixed topology or not, or whether the tri- 
angulations have to specify the piecewise linear manifold 
conditions or not. All these discussions can be summa- 
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rized in a choice of measure on the configuration space 
C . Regge calculus on a fixed triangulation has then only 
non-vanishing measure on one triangulation. In contrast 
the dynamical triangulation measure 20 is concentrated 
on triangulations with fixed edge length, that is all the 
geometric degrees of freedom are put into the choice of 
triangulation. Group field theories 21 - generate sums over 
topologies, 'generalized' triangulations and variables at- 
tached to the simplices of these triangulations. 

With the enlargement of the configuration space the 
question arises whether some form of diffeomorphism or 
more general gauge symmetries act on this space involv- 
ing configurations based on different triangulations. For 
instance the author o 22 ' 23 claim that this is actually the 
case with Romer et al.— arguing for a gauge fixing pro- 
cedure on this enlarged configuration space. 

Indeed for every solution of the Regge equations (with 
vanishing cosmological constant) one can introduce triv- 
ial refinements of the underlying triangulation and in this 
way obtain again solutions of the Regge equations. These 
refinements are trivial since the new edges or triangles do 
not carry additional curvature, hence the actual geome- 
try is invariant under these refinements. 

To be more precise, given a solution of the (four- 
dimensional) Regge equations choose a 4-simplex a and 
subdivide this simplex by placing a vertex v into this sim- 
plex. Connect the vertex v to the vertices of a. This in- 
troduces five new edges, ten new triangles and five new 4- 
simplices. To the new edges there exist a four-parameter 
space of length assignments such that the deficit angles 
at the new triangles vanish, corresponding to the possible 
positions of the new vertex v in the simplex a. Also it is 
easy to see that the Regge equations of motion 

^e t cota'=0 (19) 

are still satisfied: For the new edges the sum in (fT9")) 
involves only deficit angles of new triangles. As these are 
vanishing the equations of motion are satisfied for the 
new edges. The equation of motions for the other edges 
do not change either, again because the deficit angles at 
the new triangles are vanishing. 

Hence if one attempts to include all possible triangula- 
tions into the path integral one has to face these redun- 
dant configurations. One could try to exclude triangula- 
tions with trivial subdivisions from the path integral, for 
instance by summing only over some class of 'simplest' 
bulk triangulations that interpolate a given boundary 
triangulation 24 . Indeed this works in the (topological) 
3d caseii. To this end, however, it would be valuable to 
know more about the uniqueness properties of the Regge 
solutions: For instance, given a solution with a subdi- 
vided 4-simplex, is this 4-simplex necessarily flat? If 
this is not the case, restricting to the 'simplest triangula- 
tions' could miss out on physical degrees of freedom. So 
far some available result a 25 ' 26 i 27 indicate that solutions to 
the Regge equations are not necessarily unique and can 
vanish or appear under subdivisions or coarse graining. 



A more systematic study of the behavior of Regge solu- 
tions under subdivisions would be valuable, as this could 
be helpful for obtaining the continuum limit. If the con- 
jecture holds, that subdivisions beyond some scale (deter- 
mined by either the curvature scale of the boundary tri- 
angulation or the cosmological constant) are trivial, that 
is only add gauge equivalent solutions, it might not be 
necessary to include all triangulations into the path inte- 
gral to capture essential physical properties. It might be 
even sufficient to restrict to a fixed triangulation, which 
is fine enough, i.e. solutions should have mostly vanish- 
ing deficit angles. Summing over all kinds of length as- 
signments on this triangulation would include configura- 
tions that are equivalent to geometries defined on coarse 
grained triangulations. These terms should also arise if 
one chooses to start with another equally fine triangula- 
tion. Hence one would expect physical observables only 
including scales above the discretization scale to be inde- 
pendent from the concrete choice of this triangulation. 



IV. CANONICAL THEORIES: SPACE-TIME 
SPLIT 

There are several approaches to obtain a canonical or 
(3+1) version of Regge calculus, on the one hand in order 
to develop a numerical scheme for solving the initial value 
problem and on the other for canonical quantization. The 
main difficulty in developing a canonical scheme is firstly 
the special status of the time variable and furthermore 
the question of the existence or non-existence of gauge 
symmetries. Both points reappear in some form in loop 
quantum gravity and spin foam models, so it might be 
worthwhile to discuss these. 

There are several options to define a canonical ver- 
sion of Regge calculus. One is to start from the dis- 
crete 4-dimensional Regge action and to come up with 
a canonical, that is initial value formulation that repro- 
duces the 4-dimensional Regge solutions. This can be 
summarized as first discretizing and then performing a 
(3 + 1) -splitting. Another option is to perform first the 
(3+1) splitting and then to discretize. In the latter 
case the additional question arises whether to aim for a 
continuous-time formulation or an evolution in discrete 
time steps. 

I n 28 i 29 the second viewpoint is adopted, that is to con- 
sider the time evolution of a simplicial 3-manifold in con- 
tinuous time. Basically the (3 + 1) form of the action is 
evaluated on a three-dimensional triangulation. The re- 
sulting canonical variables are the length of the spatial 
edges, conjugated momenta (which in the version 2 ^ are 
in a non-local relation with respect to the extrinsic cur- 
vature leading to a non-local expression for the Hamil- 
tonian) and lapse and shift functions associated to the 
tetrahedra. Accordingly for each tetrahedron there is one 
Hamiltonian and three momentum constraints. Whereas 
in the continuum these constraints are first class, indi- 
cating the gauge freedom to choose lapse and shift freely, 
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this is not the case for the discretized constraints, which 
are second class. 

To obtain a consistent evolution in which the con- 
straints are preserved during time evolution one has to 
fix-2£ the values of lapse and shift, hence there is no gauge 
freedom left. 

That first class constraints turn under discretization 
into second class ones is a general phenomenon and ap- 
pears not only in Regge calculus. For instance this hap- 
pens also for the spatial diffeomorphism constraints^ in 
a lattice version of loop quantum gravity. 

As mentioned in the introduction it would be valuable 
to have a canonical formalism that matches the amount 
of symmetries inherent in the action. A strategy to en- 
sure this is to construct a canonical formalism that re- 
produces exactly the solutions of Regge calculus. To this 
end we consider 4d triangulations and as in a continuum 
have to choose a slicing of these triangulations into con- 
stant time surfaces. Such slices can be defined so that 
constant time surfaces are built from the 3d boundaries 
of the 4-simplices, and such that every vertex of the 4d 
triangulation is part of at least one constant time sur- 
face. Note that there might be many choices of possible 
slicings for a given triangulation. If the canonical for- 
malism could reproduce all such slicings it would ensure 
a discrete version of the multi-fingered time picture^! in 
continuum general relativity. 

But generically one will encounter the problem that 
constant time surfaces in the same triangulation might 
have a different number of edges and hence would lead 
to phase spaces of different size. The time evolution map 
would then have to act between different (sized) phase 
spaces. It would be of interest to have a formalism that 
can accommodate such a feature in a canonical formal- 
ism, as it might be relevant in particular for cosmological 
application. 

For the time being we restrict ourselves to cases where 
the 3d triangulation is the same for every time step. This 
restriction could lead to very rigid time evolution. Fortu- 
nately this is not the case, as in the so-called puff pastry 
evolution scheme^ an arbitrary 3d triangulation can be 
evolved in many different ways, so that consecutive 3d 
triangulations do not change. 

That is assume that we have a 4d solution of the Regge 
equations with a 3d boundary representing the equal time 
surface. This solution can be evolved further by choosing 
a vertex v in the boundary to which an edge (the 'tent 
pole') is glued. The other vertex of this edge has then 
to be connected to the vertices in the 3d boundary ad- 
jacent to v. In this way a tent is erected and we pushed 
the hypersurface locally forward in time. Since we can 
choose (apart from restrictions arising from causality in 
Lorentzian space-times) at which vertices to evolve the 
hypersurface this evolution scheme is very reminiscent of 
a discrete version of multi-fingered time. The question 
is whether we can find also a continuous or approximate 
version, that is the freedom to choose lapse and shift at 
each vertex. Indeed there are numerical indications that 



this is the cas o ' 18 . 

Evolving a vertex amounts to solving the Regge equa- 
tions associated to the new inner edges, that is the pole 
edge and the edges that are adjacent to the vertex v in the 
initial hypersurface. These equations have to be solved 
for the length of the pole edge and the length of the 
new boundary edges, hence we have as many equations 
as variables to solve for. The approximate freedom in 
choosing lapse and shift means that these equations are 
'approximately' dependent, i.e. values that satisfy all of 
these but four, will satisfy approximately the remaining 
four. 

Note that these equations for the new inner edges de- 
pend on the initial data, that is the edge length in the 
initial hypersurface and to the past of it. As in the contin- 
uum one would expect that these data satisfy constraints 
(at least in some approximate sense), and only if the ini- 
tial constraints satisfy these constraints the equations for 
the new edge length will be approximately dependent. 
These constraints describe 'initial data' on 3d hypersur- 
faces that can arise as 4d boundaries of Regge solutions. 

In the continuum the Hamiltonian and diffeomorphism 
constraints are given as the contraction of the Einstein 
tensor with the normal to the hypersurface as G a bn a n h 
and G a bn a respectively. The Regge equations associated 
to an edge e a can be understood as the double projection 
of the Einstein tensor G a be a e b to this edge. Therefore we 
will expect that one of the four constraints is given by 
the Regge equation associated to the tent pole. 

To make this more explicit we have to first find a suit- 
able set of initial data, that is phase space variables that 
can be associated to a 3d triangulation and then have to 
be evolved in discrete time steps to reproduce a Regge so- 
lution. To this end we can adopt the formalism of consis- 
tent discretization^ - as we will see it amounts essentially 
to a rewriting of the Regge equations of motion as canon- 
ical transformations on a phase space with length and 
conjugated momenta. The Regge action (with boundary 
terms) for the 4d tent region serves as a generating func- 
tion for the canonical transformation implementing the 
discrete time evolution. This will ensure that the sym- 
metries of the Regge action are reflected in the canonical 
formalism. 

The consistent discretization scheme has been applied 
to 3d Regge calculus (with non-vanishing cosmological 
constant) by Gambini and Pullin^-. But there a hyper- 
cubical lattice was chosen making the analysis very com- 
plicated. The advantage of the puff-pastry evolution 
scheme is that we can reduce our considerations to a 
very small region, in this way we have to deal only with 
a very small number of equations. A detailed treatment 
of the 4d dimensional case will appear elsewhere^, here 
we will consider a simplified 3d example corresponding 
to our earlier discussion in section [TTJ As we will see the 
exact and approximate symmetries discussed there are 
reflected in the canonical formalism. 
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A. A mini— triangulation 

The local evolution scheme describes above allows us 
to resort to a very simple example in order to illustrate 
the mechanism and the main points. Consider 3d Regge 
calculus with and without cosmological constant. (The 
results do not change for the analogous example in 4d). 
The equal time surfaces in our example will be given by 
boundaries of tctrahedra, moreover to simplify calcula- 
tions we require that the three edges of the 'base' of the 
tetrahedron have all length a and the three edges meet- 
ing at the 'tip' of the tetrahedron have all length b. We 
evolve these surfaces by the evolution described above, 
that is we join an edge (the pole) with length t to the 
tip of the tetrahedron and connect the other end of this 
edge to the vertices of the base, see picture [T] 




FIG. 1: A three-valent vertex is evolved in discrete steps by 
erecting a pole and connecting the other vertices to the end 
of the pole. The different types of line indicate the different 
time steps. 

Hence we have two configuration variables b and t, the 
edge length a does not change in time and therefore is 
merely a parameter of the system. Performing a step 
from time n to time n + 1 three simplices are added that 
hinge together at the pole. The action (including bound- 
ary terms) for this simplicial complex is given by 

S n (,b n ,b n+1 ,t n ) = -b n (tt - 26(b ni b n+ll t n )) 

-b n+ i {n - 29(b n+ll b ni t n )) 
—t n e(6„+i, b n , t n ) 
+XV(b n ,b n+u t n ) (20) 

where 9 are the dihedral angles at the spatial edges (note 
that due to the symmetry of our problem we can use 
the same functions at two consecutive time steps, but 
we have to exchange the arguments) and e is the deficit 
angle at the pole. The volume of the three tetrahedra is 
denoted by V(b n , 6 n +i, t n ). These geometric quantities 
are given by 



6(b,b',t) 
e(b,b',t) 
V(b,b',t) 



■ 



-b ,2 +t 2 ) 



2tt— 3 arccos 



V4b' 2 -aV 2 (b 2 b' 2 +b 2 t 2 +b' 2 t 2 )-b4-b'4-t4 
'2(b 2 b' 2 + b 2 t 2 + b' 2 t 2 -a 2 t 2 )-b 4 -b' 4 -^ 



2(b 2 6' 2 + 6 2 t 2 + b /2 t 2 ) _ i) 4_ i ,/4 _ 

\yf2{b 2 b' 2 +b 2 t 2 +b' 2 t 2 )-b i -b' i -t i -a 2 t 2 . (21 



We expect the length t of the time like edge to act as 
the lapse, the shift is vanishing due to our symmetry 
requirements. 

Following the consistent discretization scheme^ we 
rewrite the equations of motion for the spatial edge 
length b n and the tent pole t n 



d 

= -777— (S n -i{bn-i, b n , t n _i) + S n (b n , b n +i, t n )) 
ob n 



d 

= Tjj—S n (b n ,b n +i,t n ) 

as a canonical transformation given implicitly by 
dS(b n ,b n+ i 

1 "Tl ) 



(22) 



Pn 



Pn+1 



Pn+1 



db n 

3(ir-28(bn,b n+1 ,t n ))-\ 

dS(b n , b n +i,t n ) 
dt n 

e(bn,b n +i, t n ) — A 
dbn+i 

-3(7r-20(6 n+1 , &„,*„)) + A 

dS(b n ,b n+1 ,t n ) 
dt n +i 



(23) 



dV(bn,b n +i,t n ) 



db n 



(24) 



dV(b n ,b n+ i,t n ) 
dt n 



(25) 



dV(b n ,b n+1 ,t n ) 



dbnA 



(26) 



= 



Derivatives of the angles vanish due to the Schlafli iden- 
tity. 

Equation (|26p demanding vanishing of the t— 
momentum, is a constraint. It corresponds to the contin- 
uum result that the momentum conjugated to the lapse 
is vanishing. As in the continuum requiring that this 
constraint is preserved in time results in a secondary con- 
straint, namely that the right hand side of (|24|) should be 
equal to zero. In the continuum this secondary constraint 
is the Hamiltonian constraint. 

To obtain a constraint on the canonical data with the 
same evolution label n, we have to solve equation (|23p 
for b n+ i as a function of b n ,p n and t n and to use this 
expression in 



C x (b,p b ,t) = C o (b,p b ) + XF x (b,p\t) = . (27) 
In the case A = this constraint does not depend on t 



C (b, P b 



2tt — 3 arccos 



'cos(^)(46 2 -a 2 )-a^ 



(28) 

That the constraint depends only on b and p b expresses 
the dynamical fact that the equal time surfaces always 
bound a flat tetrahedron. This tetrahedron might be 
subdivided into smaller ones, but since the deficit angles 



8 



have to vanish the details of the subdivision do not mat- 
ter. The momenta conjugated to the variable b are given 
by the dihedral angles of this tetrahedron, which can be 
computed knowing only the edge length b (and the edge 
length a). 

The constraint (|27[) is preserved under time evo- 
lution, that is if C(b n ,Pn) vanishes this holds also 
for C(b n +i{b n ,p b n ,t n ),p b n+1 {b n ,p b n ,t n )) with t n arbitrary. 
We obtain a consistent evolution, where p b n is constrained 
as a function of b n but the lapse t n can be chosen arbi- 
trarily at every time step. 

As in the continuum where the lapse function can be 
dropped from phase space we can omit the variable t and 
are left with the edge length b and conjugated momen- 
tum p b . There is one constraint between the variables 
generating a one-dimensional gauge orbit. Hence the 
number of physical degrees of freedom is zero. The same 
result holds for arbitrary complicated triangulations (of 
spheres) in 3d gravity without cosmological constant^ 
and can also be extended to four dimensions for a special 
class of hypersurface triangulations leading to a dynam- 
ics of flat 4-dimensional spaced, see the discussion in the 
next section. 

The time steps t n can be chosen arbitrarily, that is also 
arbitrarily small. In this way we can recover a continuous 
time formalism in which the constraint ([28)1 does gener- 
ate time evolution with arbitrary choice of lapse. This 
reflects the continuous group of symmetries acting on the 
space of solutions of the covariant theory. 

For non-vanishing cosmological constant the con- 
straint (f27|) does depend on the variable t n . An analytical 
solution is not possible for A ^ 0, but one can make a 
perturbative ansatz b n+ i = ^b n+ i + A^&n+i + . . . for 
the solution of equation ([23]) and use it to obtain the 
constraint (|27p . The constraint depends on t n starting 
with the quadratic order in A. Given some initial data 
bo and Pq one can solve the constraint for to and evolve 
to the next time step where again one has to solve the 
constraint for t\ and so on. The time intervals t n are 
therefore determined by the chosen initial data bo and 
Pq. This allows to define a reduced system by solving 
the constraints at all time steps for the time intervals 
t n as functions of b n and p b n and to insert these solu- 
tions in the remaining time evolution equations. This 
has the advantage of getting rid of the constraints as well 
as some of the technical and conceptual difficulties con- 
nected to these. This 'consistent discretization' approach 
has been proposed and developed by Gambini,Pullin and 
collaborators^. 

In particular the application of this scheme to 3d Regge 
calculus (with cosmological constant) has been outlined 
ir*22, choosing a rigid slicing based on a hyper-cubical 
lattice. Note however that this strategy does not work 
in instances where the discrete action exhibits symme- 
tries, as in the case A = (even partly in four dimen- 
sions). For 3d gravity the question arises how to obtain 
the continuum limit, which has a finite number or even 
zero physical degrees of freedom, starting with a discrete 



theory with many physical degrees of freedom. Moreover 
this limit would involve an increasing number of edges 
and in the scenario without constraints therefore degrees 
of freedom. 

On the other hand from the discussion of approxi- 
mate symmetries one would expect that these symmetries 
should become exact in the continuum limit, with the 
constraints reappearing in some way. This indicates that 
choosing the (spatial) edge lengths and conjugated mo- 
menta as free variables by solving for the time like edge 
lengths may lead to a bad parametrization of the dynam- 
ically allowed configurations. The constraints (|27| could 
either have no solutions t n at all or lead to rather large 
edge lengths t n , spoiling the continuum limit. There is 
only a small range of variables (b, p b ) that will cause small 
time intervals. 

For this reason we advocate to keep the constraints 
explicitly and use them as a controlling device for the 
continuum limit. That is at the initial time step we keep 
to = e as a free parameter and find initial values bo and p b 
satisfying the constraint. These data have to be evolved 
using the equations (|23H26[) . The difference to an exact 
reparametrization invariant theory (as the A = case) is 
that now b n ,p b n and t n are fully determined by the evolu- 
tion equations for n > 0. In the case that t n ~ e, that is if 
the time intervals remain small, a (time) continuum limit 
e — > can be obtained. Numerical evolution scheme o 17 i 18 
for Regge equations follow an opposite strategy in choos- 
ing the equivalents of lapse and shift freely and ignoring 
in every time step four equations of motion which are 
therefore treated as constraints. One hopes that the vio- 
lation of these constraints does not grow during evolution 
and tends to zero in the continuum limit. Indeed numer- 
ical evidence for this scenario has been found. In this 
example the constraints are implemented exactly and the 
question is how the time steps grow. The same numerical 
evidence can be taken as an indication that if the initial 
data are chosen carefully, these remain bounded through 
evolution, although this deserves further study. 

To keep the constraints explicitly is closer to the con- 
tinuum time canonical formalisms with dynamically de- 
termined lapse and shift mentioned before and to the 
viewpoint developed in the work>2£ which argues that dis- 
cretizations will in general induce complicated gauge fix- 
ings. But the examples considered in this paper differ 
from the approach followed here. The work2& starts di- 
rectly with a Hamiltonian formalism in which there are 
explicit constraints apriori not involving the time inter- 
vals t n - Evolution of these constraints under a finite 
canonical transformation leads however to the appear- 
ance of t n in the time evolved constraints. Conservation 
of the constraints in time require the t n to be (gauge) 
fixed. Here we start from a discrete action and end up 
with constraints C(t n ,b n ,p n ) depending on the time in- 
tervals directly. 

We want to stress that this is the main difference to 
a continuous time canonical formalism, in which one as- 
sumes exact constraints, independent of lapse and shift. 
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The discrete time scenario makes the notion of 'approx- 
imate constraints' precise: we can see these as a fam- 
ily of exact constraints labeled by Lagrange multipliers. 
Choosing these Lagrange multipliers in a small interval 
(around zero) we obtain as allowed configurations for the 
remaining phase space variables a 'thickened constraint 
hypersurface'. If in the continuum limit the approximate 
symmetries are restored also the constraints should be- 
come exact: One can keep either the range for the La- 
grange multipliers fixed, then the 'thickness' of the con- 
straint hypersurface should decrease. Or we keep the 
'thickness' constant, then the range of Lagrange multi- 
pliers should increase until the full freedom of their choice 
is restored. 

In the quantum theory it is therefore consistent to im- 
pose the constraints only in an approximate sense. For 
instance in a path integral the integration over the La- 
grange multipliers should result in an exact projector on 
the constraints only if there is a solution (saddle point 
in a semi-classical approximation) to every choice of La- 
grange multiplier. A path integral with approximate 
symmetries can only impose constraints approximately 
on the canonical states. 



V. DISCRETE HYPERSURFACE 
DEFORMATION ALGEBRAS 

Dirac's hypersurface deformation algebra 

{H[M],H[N]} = D[g ab (Md b N - Nd b M)} 
{D[M a ],H[N}} = H[M b d b N] 
{D[M a }, D[N b }} = D[M b d b N a - N b d b M a ] (29) 

can be understood as the canonical incarnation^ of 
the space-time diffeomorphism group. Here H[N] — 
J s N(a)H(a)da is the Hamiltonian constraint integrated 
with a lapse function N(a) over a spatial hypersurface S 
whereas D[N a ] = J„ N a H a (a)da are the spatial diffeo- 
morphism or vector constraints integrated with the shift 
vector field N a (a). The algebra is very complicated as 
we have in the first relation of (f2"9"|) a structure function 
appearing that involves the inverse g ab of the spatial met- 
ric. 

The hypersurface deformation algebra is describing a 
universal^! relationship between infinitesimal deforma- 
tions of a hypersurface in the embedding space. The 
Hamiltonian constraints generate deformations of the hy- 
persurface normal to itself whereas the vector constraints 
generate tangential deformations, that is spatial diffeo- 
morphisms. The universality, i.e. independence from 
the field content of the theory, follows from the fact 
that the algebra can be derived by purely geometric 
considerations^!. We want to stress that therefore such 
a notion of hypersurface deformation algebra should also 
exist for discretized geometries. Also the structure con- 
stants and structure functions should be derivable by ge- 
ometrical considerations. Another question is whether it 



has universal properties and whether there are classical 
and quantum representations of this algebra realized by 
phase space generators and operators respectively. 

It has been shown^, that under certain assumptions 
including the nature and number of phase space vari- 
ables, canonical general relativity provides the unique 
phase space representation of the algebra (|29|) . A quan- 
tum representation of this algebra could therefore ensure 
the correct classical limit of the theory. Because of the in- 
tricacies and unusual form of the algebra this however has 
not been achieved so far for any theories, even free para- 
metrized field theories- 8 in space-time dimension larger 
than two. 

Given the significance of Dirac's algebra for the contin- 
uous theory it would be valuable to have a notion of a hy- 
persurface deformation algebra for discrete geometries^. 
Also here two different notions are possible, one with 
continuous and one with discrete time or deformation 
parameters. The next question is to construct (at least 
classical) representations, that is canonical theories ex- 
hibiting these symmetry algebras. Of course it would be 
interesting to know in which sense the universality fea- 
tures and the uniqueness proof can be carried over to the 
discrete cases. 

For continuous deformation parameters the easiest case 
to consider are hypersurface deformation algebras de- 
scribing embeddings of piecewise linear hypersurfaces 
into flat space. As in the continuum the algebras itself 
can be defined in any dimension by purely geometric con- 
siderations, that is by deriving the commutation relations 
of infinitesimal deformations of such hypersurfaces. One 
question would be how to define such deformations. The 
most elementary notion (and the one matching the sym- 
metries encountered so far in covariant Regge calculus) 
is to use vertex displacements, i.e. to change the position 
of one vertex in the embedding space and to accordingly 
adjust the length of the adjacent edges. Displacements 
of higher dimensional simplices can be combined out of 
vertex displacements. 

In the continuum Dirac's algebra is defined by using 
two types of geometric deformations, the (normalized) 
deformations normal to the hypersurface and tangential 
ones. Such a splitting is ambiguous at a vertex in a piece- 
wise linear geometry, so some choice for a normal to the 
vertex has to be done. The final hypersurface deforma- 
tion algebra will depend on this choice, details will appear 

Note that even a tangential vertex deformation will 
change the geometry of the spatial hypersurface (with 
the exception of intrinsically flat vertices), whereas in the 
continuum the spatial diffeomorphism constraints do not 
change the intrinsic geometry. Nevertheless for instance 
for three-dimensional gravity all three vertex deforma- 
tions are realized as constraints and are needed to obtain 
the correct counting of physical degrees of freedom for 
3d gravity. The naive viewpoint that the spatial trian- 
gulation represents the geometry of the hypersurface, i.e. 
that the spatial diffeomorphism have been already fac- 
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tored out completely, leads to an incorrect result, see 
also the discussion ini£. 

This leads us to the next question, namely if there 
exist any representations of these algebras. In three 
dimension such representations are provided by several 
canonical formulations of discretized three-dimensional 
gravit y 34 ' 41 i 42 , which can be understood as the dynamics 
of hypersurfaces embedded into flat 3d space. One for- 
mulation for (Euclidean) 3d gravity is in terms of SU (2) 
connection variables, see for instance 4 ^. The phase space 
variables are given by SU(2) holonomies and triad vari- 
ables associated to the edges of the triangulations. The 
dynamics is given by two sets of first class constraints, the 
Gaufi constraints and the so called flatness constraints. 
The Gaufi constraints can be used to reduce the phase 
space to SU(2) invariant variables giving the length of 
the edges and conjugated momenta providing a canoni- 
cal formulation of 3d Regge calculu s 34 ! 41 . (The example 
in section TlV Al with vanishing cosmological constant is a 
symmetry reduction of these theories.) The flatness con- 
straints generate vertex displacements. A priori the flat- 
ness constraints are Abelian and do not have a geomet- 
ric invariant meaning as normal or tangential generators. 
Such a split can however be introduced by using appro- 
priate linear combinations (with phase space dependent 
coefficients) of the flatness constraints. The new algebra 
is still first class and reflects the geometric commutation 
relations of the infinitesimal deformations. 

In four dimension a representation for such an alge- 
bra can be defined for a restricted class of spatial three- 
dimensional triangulations 35 . This subclass can be un- 
derstood as a sub-sector of Regge calculus, whose trian- 
gula tions only allow flat solutions (see also the papers^ 3 - 
for defining a flat geometric sub-sector in a canonical 
formulation of BF-theory) . The reason for the restricted 
class of triangulations is the following: A 2d triangulation 
with prescribed edge lengths can be generically (ignoring 
triangle inequalities, self intersections and similar issues) 
embedded into 3d flat space. This reflects the dynamical 
content of 3d gravity: Given a 2d metric the constraints 
can be solved (locally in phase space) for the momenta 
which define an embedding of the corresponding surface 
in flat 3d space. 

A generic 3d triangulation with prescribed length can- 
not be embedded into flat 4d space, as this would require 
constraints on the length variables 35 . However there is 
a special class of triangulations, such as the boundary of 
a 4-simplex, which can be embedded into 4d flat space. 
More complicated triangulations with this property can 
be generated by applying (1 — 4) Pachner moves to this 
triangulation, i.e. by placing a vertex into a tetrahedron, 
subdividing it thus into four new tetrahedra. A more 
complete discussion can be found in the paper— here we 
will just give the canonical theory describing the 'dynam- 
ics' of a 4-simplex boundary embedded in flat space. 

The canonical variables are given by the ten edge 
length l e denoting vertices and conjugated momenta p e . 
A canonical transformation can be performed (in local 



patches of phase space) to area variables At (I) and con- 
jugated momenta pt = 9 g^ Pe- The dynamics is 
then defined by the ten constraints 

C t :=p t -e t (l) (30) 

where 9t is defined to be the dihedral angle at the trian- 
gle t in a flat 4-simplex. Due to the Schlacfli identity the 
constraints are first class and hence generate gauge trans- 
formations. These gauge transformations correspond to 
the vertex deformations of the 4-simplex. Since there are 
5 vertices in a 4-simplex we have 20 possible deforma- 
tions, however there are 10 combinations corresponding 
to rigid translations and rotations of the 4-simplex which 
do not change the geometry of the 4-simplex. 

Taking the following linear combinations of constraints 

^£^.--£^> <*> 

we obtain generators that change only one length of one 
edge at a time (and therefore commute). The constraint 
C e induce therefore a displacement of either of the ver- 
tices in the direction normal to the other edges adjacent 
to that vertex. (The ambiguity in associating a vertex 
stems from the fact that a transformation at both ver- 
tices amounts to a rigid transformation of the 4-simplex) . 
Taking linear combinations of these constraints we can 
define generators leading to arbitrary vertex displace- 
ments and compute a deformation algebra. 

Note that in equation (|3"D|) the triangle areas seem to 
appear as an easier set of variables to deal with. The 
conjugated momenta are directly related to the dihedral 
angles, encoding the extrinsic curvature of the triangu- 
lated 3d surface. Indeed also in loop quantum gravity ar- 
eas and 3d dihedral angles between triangles arise as the 
simplest kind of geometrical variables. These variables 
provide however an over-complete set and need therefore 
to be constrained in order to yield consistent length as- 
signments to the edges of the triangulations. Fortunately 
such constraints are available for 4d triangulations 4 ^ and 
can be used to define an area-angle Regge action. The 
same constraints arise for 3d triangulations 3 — and can 
therefore be used in a canonical analysis of the area-angle 
Regge action 4 ^. The advantage of such a formalism is its 
closeness to spin foams and loop quantum gravity. In fact 
the variables used there are related by a Gaufi constraint 
reduction 35 to the area and angle variables. 

The incidences where one can define representations of 
an (infinitesimal) hypersurface deformation algebra cor- 
respond exactly to the known appearances of non-trivial 
symmetries in the corresponding covariant versions, that 
is Regge calculus. If there are no corresponding symme- 
tries in the action, it is questionable that a first class con- 
straint algebra reflecting infinitesimal hypersurface de- 
formations of discretized hypersurfaces can be defined 
exactly, as this would generate solutions not present in 
the covariant theory. It might exist in an approximate 
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sense, that is reflecting the algebra ([29]) modulo anoma- 
lous terms vanishing in the continuum limit. This has 
been shown for the spatial diffeomorphism constraints in 
a lattice approach to loop quantum gravity 30 . 

In a quantum theory a deformation algebra might 
rather involve finite transformations. Since the spectra 
of spatial geometric quantities in loop quantum gravity 
are discrete^ 5 - there might appear a similar discreteness 
in time direction. In the examples discussed so far the 
choice at which vertices to apply a tent move in the puff- 
pastry evolution scheme^ can be understood as choosing 
lapse to be either equal to zero or equal to some dynam- 
ically determined value. In a quantization such a move 
could give rise to a finite (unitary) evolution instead of a 
Hamiltonian viewed as infinitesimal generator. Similarly 
the spatial diffeomorphism constraints in loop quantum 
gravity do not exist as infinitesimal generators but only 
as finite unitary operators. 

This understanding of the hypersurface deformation 
algebra corresponds to different choices of triangulations 
in the covariant theory. The puff pastry evolution scheme 
will generate a certain subclass of triangulations. A com- 
plete equivalence can be obtained if one considers all the 
allowed Pachner moves in the (n — l)-dimensional hy- 
persurface. These can be related to gluing or removing 
n-simplices from the hypersurface thus building up an n- 
dimensional triangulatio n 34 i 46 . The discrete time notion 
of hypersurface deformation algebras therefore depends 
on a better understanding of any symmetries associated 
with the choice of triangulations. 



VI. LOOP QUANTUM GRAVITY AND SPIN 
FOAMS 

In this section we will consider implications of the re- 
sults collected in the previous section for quantum grav- 
ity models, in particular spin foams and canonical loop 
quantum gravity 5 -*^. 

Similar to quantum Regge calculus spin foam models 
are based on a path integral over discretized geometries. 
What differs from (traditional) quantum Regge calcu- 
lus is the choice of variables and the discrete integra- 
tion measure over these variables . This leads to a sum 
over spin labels (ST/(2)group representations) and group 
intertwiners, associated to subsimplices of the triangu- 
lation. The discrete integration measure is motivated 
among other things by canonical loop quantum gravity, 
where the spectrum of kinematical geometric operators 
is discrete and encoded in such spin labels and intertwin- 
ers. To obtain the amplitude associated to a given set of 
labels one considers the (discretized) Plebanski action^ 
that formulates general relativity as a second class con- 
strained BF theory. The constraints are imposed by dif- 
ferent methods 4 ^ onto the quantum amplitudes. Further- 
more the amplitudes depend on a choice of discretized 
measure which is often left open to some extend in the 
definition of spin foam models. 



So far the strategy to obtain a semi-classical descrip- 
tion is to consider a fixed triangulation and to consider 
the large spin limit. This is a limit where the discretiza- 
tion scale given by the spectra of geometric operators, 
i.e. the Planck length becomes small compared to edge 
length. (There is a second discretization scale defined 
by the underlying piecewise linear triangulation, in other 
words the typical edge length.) In this limit one hopes 
to find an effective description given by the Regge action 
and indeed there are indications that certain models sat- 
isfy this criterion 50 . In so far one can understand spin 
foam models as a quantization of Regge calculus. For 
other models constructed in a more direct way from first 
order versions of Regge's action see the work 5 ^. The ques- 
tion arises whether the exact or approximate symmetries 
of the discretized action are reflected in spin foam mod- 
els, that is in the amplitudes defined by the action and 
the choice of measure. 

There is another viewpoin t 52 ' 53 of the path integral 
which makes it absolutely necessary to study the symme- 
tries of the path integral: A path integral for constrained 
theories such as gravity can be understood as a projec- 
tor onto physical states, i.e. wave functions that satisfy 
the Hamiltonian and diffeomorphism constraints. (This 
is not a proper projector as it will map normalizable wave 
functions to unnormalizable ones.) 

The heuristic idea behind is that the path integral im- 
plements an averaging over the gauge orbits generated 
by the constraints, resulting in wave functions that are 
constant along these orbits. This can only work, if the 
path integral respects the symmetries, in particular if 
the measure is gauge invariant. So far this idea is rather 
heuristic for 4d gravity and has been made only rigor- 
ous for 3d gravity 5 - 3 -, where the path integral has all the 
necessary symmetries 1 ^. 

As was mentioned before a discretized path integral 
might have a complicated symmetry structure: diffeo- 
morphism or translation symmetries associated to trivial 
vertices can be implemented exactly, for non-trivial ver- 
tices there might be only approximate symmetries. Here 
it would be valuable to know whether this can lead at 
least to an approximate 'projector', that is a Gaussian 
compared to a delta function. This would be consistent 
to the discussion in section [TV Ai where we pointed out 
that in a discrete theory the allowed initial data are de- 
scribed by a 'thickened' constraint hypersurface. In the 
continuum limit, in which the symmetries become exact, 
the projector should also become exact. Approximate 
symmetries are reflected by very small eigenvalues in the 
Hessian of the action, which should converge to zero in 
the continuum limit. The integral over the directions 
corresponding to these eigenvalues is an integral over the 
Lagrange multipliers parameterizing the (approximate) 
gauge orbit. If the action and the path integral mea- 
sure is approximately constant along these gauge orbits 
it should result in a projector onto approximately gauge 
invariant states. This behavior should be studied in toy 
models. Other ideas how to construct a projector on 
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physical states will be mentioned below. 

So far we discussed the symmetries of the action, but to 
obtain a 'projector onto physical states' also the path in- 
tegral measure has to respect the symmetries. In the con- 
text of spin foam models this has been discussed exten- 
sively by Bojowald and Perea^. They derive two kinds of 
requirements. The first requirement can be understood 
to arise from the continuum diffeomorphism symmetry 
leading to vertex translations in the discretized action. 
Together with the second class constraints appearing in 
the Plebanski action the constraints associated to the 
vertex translations (either second or first class) form a 
complicated constraint algebra, which has been so far 
only computed in the continuum case^. This constraint 
algebra is however needed to compute the Fadeev-Popov 
determinant to obtain the correct path integral measure. 

The second requirement is related to symmetries which 
involve a change of triangulation and labels (so called 
moves). In spin foam models labels are allowed to have 
zero spin leading to vanishing edge length - therefore 
there are more moves that leave the encoded state un- 
changed. The amplitudes associated to two configura- 
tions related by such moves should be the same. This 
leads to a number of consistency conditions for the am- 
plitudes, in particular restricting the choice of face and 
edge amplitudes, that are often left open to choice. 

Furthermore Bojowald and Perez argue that a con- 
sistent choice of amplitudes has to lead to divergen- 
cies. This has to be expected as the path integral in- 
cludes a sum over non-compact gauge orbit. For the 3d 
spin foam model this has been discussed by Freidel and 
Louapre^ and the symmetries gauge fixed. For the 4d 
case, if the spin foam model reproduces the symmetries 
of the Regge action, one would expect at least diver- 
gencies associated to the exact translation symmetry of 
trivial vertices, whereas approximate symmetries associ- 
ated to non-trivial vertices could still lead to finite results 
(however diverging in the continuum limit). Here trivial 
vertices are those for which saddle points in the Regge 
action exist such that the deficit angles at adjacent trian- 
gles vanish. This depends on the boundary conditions for 
generic vertices but will always be the case for five valcnt 
vertices. 

Spin foam models are hoped for to provide a projec- 
tor on diffeomorphism invariant states in canonical loop 
quantum gravity. On the other hand such states are de- 
fined to be annihilated by the diffeomorphism and Hamil- 
tonian constraints, and we will discuss certain aspects of 
these constraints in the remainder of this section. 

Loop Quantum gravity provides a rigorous construc- 
tion of a kinematical Hilbert space which carries a repre- 
sentation of (kinematical) geometric operators. As was 
already mentioned these turn out to have a discrete^ 
spectrum, which is one reason for the choice of a discrete 
path integral measure for the spin foam models. This 
makes it very tempting to associate a discretized piece- 
wise linear geometry to the eigenstates, the so-called spin 
networks, of these operators, a picture close to the con- 



struction of spin foam models as discretized path inte- 
grals. 

But loop quantum gravity is a quantization of the con- 
tinuum theory and does not proceed by considering dis- 
cretized classical geometries and then quantizing. Hence 
the Hilbert space spanned by spin networks is the state 
space associated to continuum geometries. As such one 
would expect the implementation of the continuum con- 
straints and continuum Dirac algebra. Indeed the cele- 
brated FLOST uniqueness theorem^, showing that loop 
quantum gravity provides a unique kinematical set-up 
satisfying certain assumptions, requires a representation 
of the continuous (more specifically semi-analytic) spa- 
tial diffeomorphisms. 

Another route is to first discretize the geome- 
try, obtain a canonical description and then quan- 
tize using loop quantum gravity techniques, see for 
instanc o 35 ' 43 i 58 ' 59 i 60 ' 61 i 62 where this approach is followed 
or advocated for. Here one would expect the implemen- 
tation of the (possibly second class) constraints for a dis- 
cretized geometry. The relations on the kinematical and 
dynamical level between these two approaches are not 
quite clear— and should be explored in more detail. 

As spin foam models are derived from a discretized 
path integral one would expect that boundary states of 
this path integral are states obtained from the second 
route of canonical quantization. Also at least before the 
continuum limit spin foams should project onto the phys- 
ical states defined by the constraints in the explicitly dis- 
cretized theory. Group field theories generate sums over 
all possible triangulations (or more general complexes), 
in a certain sense these theories should include informa- 
tion about continuum geometries and dynamics. Hence 
a canonical formulation of these could be connected to 
a continuum canonical theory related to loop quantum 
gravity. The immediate difficulties in obtaining such 
canonical formulations as well as its potential significance 
are sketched in the report^ 3 .. The full set of symmetries 
in the covariant version has still to be identified. 

For the continuum theory of loop quantum gravity 
a uv-finite quantization of the Hamiltonian constraints 
has been obtained by Thiemann^. The status of the 
quantum Dirac algebra in the (continuum) loop quantum 
gravity is however shrouded by several severe technical 
difficulties— i 66 ' 67 . That is the algebra of the Hamiltonian 
constraints is anomaly free, even Abelian. Abelianess 
does not reproduce the constraint algebra (|29|) but is con- 
sistent as the Hamiltonian constraints require for their 
continuum limit the topology provided by the Hilbert 
space of spatially diffeomorphism invariant states. On 
this space the diffeomorphism constraints (on the right 
hand side of the classical Poisson bracket between two 
Hamiltonians) vanish by definition. A further difficulty 
in obtaining an infinitesimal algebra is that the spatial 
diffeomorphism constraints are not quantized as infinites- 
imal generators but are only implemented via a family of 
non-weakly continuous unitary operators. For extended 
discussions of these aspects we refer to 65 ' 66 . 
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For the second route of first discretizing and then quan- 
tizing there is not a full classical definition of the con- 
straint algebra (apart from the sub-sector yielding flat 
dynamics) yet, but this might be obtained with the tech- 
niques described in section IIV Al As discussed there 
the results might be a complicated mix of second and 
first class constraints. Additionally one would expect 
that with decreasing discretization scale the approximate 
symmetries become exact and hence the second class con- 
straints converted into first class. Before quantization 
one should test whether this is indeed the other- 
wise 'classical anomalies' are bound to appear. 

If one follows the approach to quantize the discrete the- 
ory and then to perform the continuum limit the question 
arises of how to deal with the constraints. Usually sec- 
ond class constraints have to be taken care off at the 
classical level as they cannot be implemented exactly on 
quantum states: The uncertainty relations require the 
fluctuations of the second class constraints to be of the 
order of the expectation value of the commutator. Again 
one might hope that heuristically the expectation values 
of this commutator should become small with the dis- 
cretization scale (encoded in the state used to obtain the 
expectation value). 

To solve the second class constraints at the classical 
level and to introduce Dirac brackets would be overly 
complicated and might not allow for a continuum limit. 
Therefore methods are asked for that can deal with mixed 
constraint systems on the quantum level. There are sev- 
eral such proposal o 61 ' 67 i 68 based on different motivations. 
But they have in common to consider just one ('master') 
constraint, typically a sum or integral of the squares of 
the original constraints. This constraint is used to con- 
struct a projector on the physical Hilbert space, either 
via path integral methods^ or via spectral analysis. In 
the latter case, since the master constraint operator has 
to be positive definite, a physical Hilbert space can be 
defined as the (generalized) eigenspace corresponding to 
the minimum of the spectrum. (Alternatively one can 
consider the generalized eigenspaces from a small inter- 
val around zero and in this way obtain a proper projec- 
tor. This would correspond to the 'thickened constraint 
hypersurface' picture from section IIV AH As there are 
second class constraints this minimum will not be zero 
and may or may noi£Z- converge to zero in the contin- 
uum limit. The 'uniform discretization' method— has 
been applied to a spherically reduced model with spa- 
tial diffeomorphism constraints and it is found that the 
minimum of the spectrum converges to zerc 69 . 

Following these approaches a physical Hilbert space 
can be constructed even with second class constraints. 
The question arises how to detect quantum anomalies, 
i.e. to test whether the (hypothetical) quantized sec- 
ond class constraints become first class in the continuum 
limit. The minimum of the spectrum going to zero in the 
continuum limit might be one indication but there are 
examples^ 7 - based on non-compact Lie groups in which 
the minimum is non-zero despite a non-anomalous quan- 



tum constraint algebra, therefore one might have to allow 
for non-zero values. Furthermore in a classical reduction 
second class constraints lead to the elimination of only 
half the number of degrees of freedom as for first class 
constraints. One has to make sure that unphysical gauge 
degrees do not survive the continuum limit. 

Quantum anomalies might depend on operator order- 
ings and the precise definition of the master constraints 
but could also arise do to a physical reason, one possibil- 
ity being a fundamental discretization scale for time like 
intervals. To differentiate between these issues it might 
be necessary to quantize the original constraints and to 
study the commutator algebra on the family of physical 
Hilbert spaces which arise in taking the continuum limit. 



VII. DISCUSSION 

Incorporating diffeomorphism symmetry into quantum 
gravity models is a very difficult task. It can however 
ensure the correct semi-classical and continuum limit of 
quantum gravity. We expect that models involving a dis- 
cretization of the action will break at least a part of the 
symmetries to approximate ones, although existent dis- 
cretizations, as for instance of the Plebanski action, have 
yet to be checked. One method is to consider the eigen- 
values of the Hesse matrix evaluated at a solution. A 
minimal requirement for discretized actions is that sym- 
metries corresponding to translations of flat vertices are 
preserved exactly. 

This opens the question whether one can construct bet- 
ter actions or discretizations in which all the symmetries 
are exact. Indeed this is possible for 3d gravity with 
cosmological constant 70 . For 4d gravity renormalization 
flow methods could yield such (possibly non-local) ac- 
tions. Lehto et.alj 7 ^ consider an action with quadratic 
curvature terms in order to force the dynamics to give 
nearly flat geometries, as for flat geometries the symme- 
tries are exact. In addition symmetries corresponding 
to a free choice of triangulation could appear in covari- 
ant and in canonical formalisms. It would be of inter- 
est whether symmetries involving only a change of labels 
and symmetries involving also a change of triangulation 
or more general simplicial complex can be related. 

On the other hand one can proceed with actions where 
the symmetries are only approximate. The methods of 
consistent discretization can be used to obtain canoni- 
cal theories that reflect the exact and approximate sym- 
metries of such actions. Approximate symmetries lead 
to approximate constraints - these are families of con- 
straints labeled by Lagrange multipliers. The algebra of 
such constraints should be studied and it remains to be 
seen whether the continuum constraint algebra can be 
recovered in the limit. At least for discretized hyper- 
surfaces embedded in flat space a first class algebra of 
constraints can be obtained. 

It is necessary to have the same symmetries in the co- 
variant action and the canonical formalism in order to 
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make the idea of a path integral acting as a projector 
on physical states consistent. To this end also the path 
integral measure should be taken into account. Here con- 
ditions for an approximate invariance have to be formu- 
lated. Alternatively a projector on physical states can 
be constructed using canonical quantization and in most 
cases some kind of master constraint. Results have to 
checked carefully for anomalies. 

Boundary states^ arising from path integrals over a 
region with boundary should reflect the symmetries of 
the path integral. The same holds for observables act- 
ing on such states. As most of the symmetries are bro- 
ken this might not constrain states and observables as 
much as requiring exact gauge invariance. Note how- 
ever that non-gauge invariant variables might be very 
sensitive to discretization ambiguities and furthermore 
reflect unphysical degrees of freedom, as is emphasized 
in the consistent discretization program^. Indeed in the 
limit in which the constraints arc imposed exactly on the 



boundary states, operators that do not commute with 
the constraints will have unbounded fluctuations due to 
the uncertainty relations. Hence observables commuting 
either approximately or exactly with the constraints are 
required. As an exact commutation might restrict the 
space of observables to much for anomalous constraint 
algebras, proper conditions have to be formulated. Tech- 
niques to compute observables are available for contin- 
uum theories^ and should be adapted to discrete ones. 
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